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Abstract 

A method to determine the running of a from a measurement of small-angle Bhabha scattering 
is proposed and worked out. The method is suited to high statistics experiments at e + e~ colliders, 
which are equipped with luminometers in the appropriate angular region. A new simulation code 
predicting small-angle Bhabha scattering is also presented. 
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1 Introduction 



The electroweak Standard Model SU(2) U(l) contains Quantum Electrodynamics (QED) as a 
constitutive part. The running of the electromagnetic coupling a is determined by the theory as 

a{q) = l-Aa(^) ' (1) 

where a(0) = ao is the Sommerfeld fine structure constant, which has been measured to a precision of 
3.7xl0 -9 [Q; Aa(q 2 ) positive arises from loop contributions to the photon propagator. The numerical 
prediction of electroweak observables involves the knowledge of a(q 2 ), usually for q 2 ^ 0. For instance, 
the knowledge of a(m%) is relevant to the evaluation of quantities measured by the LEP experiments. 
This is achieved by evolving a from q 2 = up to the Z-mass scale q 2 = m 2 z . The evolution expressed 
by the quantity Aa receives contributions from leptons, hadrons and the gauge bosons. The hadronic 
contribution to the vacuum polarization, which cannot be calculated from first principles, is estimated 
with the help of a dispersion integral and evaluated [2] by using total cross section measurements of 
e + e~ — > hadrons at low energies. Therefore, any evolved value a(q 2 ) particularly for \q 2 \ > 4mJ, 
is affected by uncertainties originating from hadronic contributions. The uncertainty on a(m^)~ 1 
induced by these data is as small as ±0.09 nevertheless it turned out [3j that this limits the accurate 
prediction of electroweak quantities within the Standard Model, particularly for the prediction of the 
Higgs mass. 

While waiting for improved measurements from BEPC, VEPP-4M and DAFNE as input to the 
dispersion integral, intense efforts are made to improve on estimating the hadronic shift Aahad, as 
for instance @] - 0, and to find alternative ways of measuring a itself. Attempts have been made to 
measure a{q 2 ) directly, using e + e~data at various energies, such as measuring the ratio of e + e~7/e + e~ 
[S] or more directly the angular distribution of Bhabha scattering [3] . 

In this article the running of a is studied using small-angle Bhabha scattering. This process 
provides unique information on the QED coupling constant a at low space-like momentum transfer 
t = — \q 2 \, where 

t= -- s (l-cos0) (2) 

is related to the total invariant energy y/s and to the scattering angle 8 of the final-state electron. 
The small-angle region has the virtue of giving access to values of a(t) without being affected by 
weak contributions. The cross section can be theoretically calculated with a precision at the per mille 
level. It is dominated by the photonic t channel exchange and the non-QED contributions have been 
computed ^U] and are of the order of 10 -4 (see table QJ; m particular, contributions from boxes with 
two weak bosons are safely negligible. 

In general, the Bhabha cross section is computed (see sect. 3) from the entire set of gauge-invariant 
amplitudes in both the s and t channels. Consequently, two invariant scales s and t govern the process. 
The different amplitudes are functions of both s and t and also the QED coupling a appears as a(s) 
resp. a(t) However, the restriction of Bhabha scattering to the kinematic regime of small angles 
results in a considerable simplification, since the s channel then gives only a negligible contribution, as 
is quantitatively demonstrated in tabled Thus, the measurement of the angular distribution allows us 
indeed to verify directly the running of the coupling a(t). For the actual calculations, 9 m e / -Ekeam 
and -Ebcam 3> w e must be satisfied (see sect. 4.1). Obviously, in order to manifest the running, the 
experimental precision must be adequate. 

This idea can be realized by high-statistics experiments at e + e~ colliders equipped with finely 
segmented luminometers, in particular by the LEP experiments, given their large event samples, by 
SLC and future Linear Colliders. The relevant luminometers cover the t-range from a few GeV 2 to 
order 100 GeV 2 . 

The i-dependence of the quantity Aa(t) (eqs. ^ 121 a t small values of t is illustrated in fig. 2] It 
shows the predicted running of a in the relevant space-like region. The figure is obtained using the 
program alphaQED by Jegerlehner 0. At low energies (see fig.EJ) Aa is dominated by the contribution 
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Figure 1 : Aa versus in units of GeV in the space-like region. 



from the leptons, while with increasing energy also the contribution involving hadrons gets relevant. 
The region where hadronic corrections are critical is contained in the considered i-range. 



The experimental determination of the angular distribution of the Bhabha cross section requires the 
precise definition of a Bhabha event in the detector. The analysis follows closely the procedure adopted 
in the luminosity measurement, which is described in detail, for instance in ref. |12| . and elaborates on 
the additional aspect related to the measurement of a differential quantity. To this aim the luminosity 
detector must have a sufficiently large angular acceptance and adequate fine segmentation. The 
variable t (eq. [5J is reconstructed on an event-by-event basis. 

The method to measure the running of a exploits the fact that the cross section for the process 
e + e~ — * e + e~ can be conveniently decomposed into three factors : 



as worked out in sect. 3. All three factors are predicted to a precision of 0.1% or better. The first factor 
on the right-hand side refers to the effective Bhabha Born cross section, including soft and virtual 
photons according to ref. ^Oj, which is precisely known, and accounts for the strongest dependence 
on t. The vacuum-polarization effect in the leading photon t channel exchange is incorporated in the 
running of a and gives rise to the squared factor in eq. [3] The third factor, Ar(t), collects all the 
remaining real (in particular collinear) and virtual radiative effects not incorporated in the running 
of a. The experimental data after correction for detector effects have to be compared with eq. |3 The 
t dependence is rather steep, thus migration effects may need attention. 

This goal is achieved by using a newly developed program based on the already existing semiana- 
lytical code NLLBHA ^01 ED- A detailed description of this code called SAMBHA is given in sect. 4. 



2 The method 
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Figure 2: Contributions to Aa from leptons (dashed curve) and hadrons (solid curve) versus y— t in 
units of GeV. 



It is convenient to confront the fully corrected measured cross section with the Bhabha cross section, 
including radiative corrections in the factorized form given by eq. 01 The physical cross section is 
infrared safe JU]. This decomposition is neither unique nor dictated by a compelling physical reason; 
rather it allows the separation of the different sources of t dependence in a transparent way without 
introducing any additional theoretical uncertainty. The various factors are discussed one by one in 
the following subsections. 

3.1 The cross section da°/dt 

The differential cross section da /dt is defined as : 



The factor da B j dt is the Bhabha cross section in the improved Born approximation, which, by defini- 
tion, includes the running of a. As seen explicitly in the formulae below (eq.[SJ the term a(t)/a(0) is 
not factorized completely in the improved Born cross section. In order to have the factorized form of 
eq. |3 the t channel contribution to the running of a has been taken out. In this way, da°/dt contains 
not only the usual Born t dependence, i.e. 1/t, but also some weaker t dependences arising from s 
channel amplitudes with vacuum polarization effects taken into account |1(J| . although numerically 
small as mentioned above. 

The improved Born cross section for Bhabha scattering within the electroweak Standard Model is 
precisely known (see refs. ^JElEj)- The differential cross section da B /dt differs from da /dt by the 
inclusion of those radiative corrections that affect only the propagator of the exchanged photon. They 
form a gauge-invariant subset of all radiative corrections and are shown explicitly. It is convenient 
to decompose da B /dt into the contributions arising from the t channel (B t ), the s channel (B s ) and 
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their interference (-Bj): 



da 1 



dt 2s 



= ^He{Bt + B, + B(}, (5) 



where 
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e 2 + glf + (1 + c) 2 ( 5 " + 5a + 65 2 3 2 ) 
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X 



1 



s — m 2 + im^r^ sin 2^ 
t 1 



t — to^j sin 29 u 



g a = —z, g v = +2sin 2 e»„,), 

* = (Pi - «i) 2 = s C 1 ~ c )> 
c = cos#, 9 = PiQi- 

Here s is the total squared invariant mass, 6 W the electroweak mixing angle and the scattering angle 
between the initial and final electron with momentum p x and q 1 respectively (see ref. |10|). 

In table H the cross sections are given in nanobarns for the pure QED and electroweak cases. 
QED f denotes the contribution of the t channel pure QED Feynman diagrams. The cross sections 
are integrated over two relevant angular ranges. The table shows that the t channel photon exchange 
dominates the cross section at small angles and justifies why the process is suited for investigating the 
t dependence, and so the running of a(t). 

By comparing the values of the electroweak cross section with the pure QED one, it is seen that 
the Z-boson exchange gives a negligible contribution to small-angle scattering. In the last two lines 
(EW+VPt and EW+VP) there are numbers for the cross section with vacuum polarization taken into 
account in the t channel only, and in all channels, correspondingly. One can see that the effect of s 
channel vacuum polarization is small, as a result of the smallness of the s channel photon-exchange 
contribution itself. The last line in the table corresponds to the complete formula in eq. 

3.2 The running of a 

In eq. El the two-point functions II(£) = Aa(t) and II(s) = Aa(s) are responsible for the running of 
a in the space-like and time-like regions. In the language of Feynman diagrams the effect arises from 
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Table 1: Various cross sections in nb as a function of the centre-of- mass energy in GeV integrated 
over the two angular ranges 45-110 mrad and 5-50 mrad. The index t denotes the contribution of the 
corresponding t channel Feynman diagrams alone. The last columns are of interest for furture Linear 
Colliders. 



(GeV) 


91.187 


91.2 


189 


206 


500 


1000 


3000 








45 mrad < 9 < 


110 mi- 


ad 




y/(-t) (GeV) 


3.4 


3.4 


7.1 


7.7 


18.8 


37.5 


112.6 


QED 


51.428 


51.413 


11.971 


10.077 


1.7105 


0.42763 


0.047514 


QED t 


51.484 


51.469 


11.984 


10.088 


1.7124 


0.42809 


0.047566 


EW 


51.436 


51.413 


11.965 


10.072 


1.7105 


0.42871 


0.049507 


EW+VPi 


54.041 


54.018 


12.743 


10.745 


1.8590 


0.47303 


0.055748 


EW+VP 


54.036 


54.013 


12.742 


10.744 


1.8588 


0.47296 


0.055742 








5 mrad < 8 < 


50 mrad 




y/(-t) (GeV) 


1.1 


1.1 


2.2 


2.4 


5.8 


11.6 


34.8 


QED 


4963.4 


4962.0 


1155.4 


972.54 


165.08 


41.271 


4.5857 


QED t 


4963.5 


4962.1 


1155.4 


972.57 


165.09 


41.272 


4.5858 


EW 


4963.4 


4962.0 


1155.4 


972.53 


165.08 


41.272 


4.5885 


EW+VP t 


5075.0 


5073.5 


1190.6 


1003.3 


172.51 


43.647 


4.9603 


EW+VP 


5075.0 


5073.5 


1190.6 


1003.3 


172.51 


43.646 


4.9605 



fermion-loop insertions into the virtual photon lines: 

n(t) = + 5 



7T V 3 9 



lv)(i i+ « 3 »-| 



3 /™2 



where 



o 2 

ln^, Q 2 = -i, C(3)= 1.202 
mi 



and where the leading part of the two-loop contribution to the polarization operator is taken into 
account. The most significant part arises from the electrons and is L/3 — 5/9. 

The O(a) and 0(a 2 ) leptonic vacuum polarization has been known for many years [17| . The third- 
order (three-loop) leptonic contributions II^ 3 ) (i) have recently been calculated ^B] ■ In the Standard 
Model, 5t contains contributions from muons, r-leptons, IV-bosons and hadrons : 

S t = 5? + 61 + 5f + Sf, 
S s = 5 t (t — > s), 

which means that 5 S is obtained from 5t by substituting s by t, see ref. |10j . The contributions from 
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the leptons (1 = fi,r) and from the W are theoretically calculable and given by: 



S\ 



Vl 



v w = 




For Q 2 ^> m 2 the formula simplifies to 

The hadronic contribution cannot be calculated theoretically; instead, it can be expressed as a 
dispersion integral involving experimentally measured e + e~ cross sections: 

For numerical calculations, hadronic contributions as included in the parametrisation of refs. |Hll9j 
are adopted. 

This procedure, as usually assumed (see e.g. 0), is based on the analyticity of the function a(q 2 ) 
in the complex plane, except possibly at the energies corresponding to the Landau pole. For the 
leptonic contributions 5^'^' T this assertion is true, while for the hadronic contribution i5j lad it relies on 
the dispersion approach to the entire, non-perturbative, hadronic physics (see eq. |HJ). This ends up in 
a single analytical function that can be used to deal with the vacuum polarization in the t channel. 

3.3 The radiative factor 1 + Ar(t) and neglected terms 

For the present investigation of the small-angle Bhabha cross section only the corrections consistently 
needed to maintain the required accuracy are kept. All these corrections are included in the new code 
SAMBHA. All the following contributions have been proved to be negligible ^U] and are dropped : 

• Any electroweak effect beyond the tree level, for instance appearing in boxes or vertices with Z° 
and W bosons, running weak coupling, etc. 

• Box diagrams at order a 2 and larger 

• Contributions of order a 2 without large logarithms, leading from order a 4 (i.e. a 4 L 4 , ■ • ■ ) and 
subleading higher order (a 3 L 2 , a 4 L 3 , • • • ) 

• Contributions from pair-produced hadrons, muons, taus and the corresponding virtual pair 
corrections to the vertices (estimated to be of the order of 0.5xl0~ 4 ). 

The radiatively corrected Bhabha cross section is denoted by da/dt. Numerically it differs from 
da B /dt by less than a few per cent for small angles, depending on energy and final-state selection 
procedure. 
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4 Monte Carlo codes and comparison 



The precise determination of the luminosity at e + e~ colliders is a crucial ingredient to obtain an 
accurate evaluation of all the physically relevant cross sections. They necessarily have to rely on 
some reference process, which is usually taken to be the small-angle Bhabha scattering. Given the 
high statistical precision provided by the LEP collider, an equally precise knowledge of the theoretical 
small-angle Bhabha cross section is mandatory. In the 1990's the substantial progress in measuring 
the luminosity reached by the LEP machine has prompted several groups to make a theoretical effort 
aiming at a 0.1% accuracy |12l I2()| . This goal has indeed been achieved by developing a dedicated 
strategy. For the first time small-angle Bhabha scattering was evaluated analytically, following a 
new calculation technique JU| that yields the required precision. Analytical calculations have been 
combined with Monte Carlo programs in order to simulate realistically the conditions of the LEP 
experiments. 

The analytical results evaluated for the various contributions to the observed Bhabha cross section 
in ref. |l(Jj . were implemented into the semi-analytical code NLLBHA (for a short write-up see in 
ref. ^31)- The important feature of this code consists in the systematic account of all QED radiative 
corrections required to reach the per mille precision. On the other hand, the simulation of realistic 
experimental acceptances can only be achieved with Monte Carlo techniques. For this purpose a 
Monte Carlo code, LABSMC, was developed [2H 123 ES] 

4.1 SAMBHA-NLLBHA 

The program LABSMC, which was intended to describe large-angle Bhabha scattering at high energies, 
has been complemented with a set of routines from NLLBHA so as to be applicable to small-angle 
Bhabha scattering. This implied the insertion of the relevant second-order next-to-leading radiative 
corrections (0(a 2 L)) in the Monte Carlo code 1 , which are crucial to achieve the per mille accuracy. 
The extension to cover small angles resulted in the new code SAMBHA containing the previously existing 
features together with the following new characteristics : 

• the complete electroweak matrix element at the Born level; 

• the complete set of 0{a) QED radiative corrections (including radiation from amplitudes with 
Z-boson exchange); 

• vacuum-polarization corrections by leptons, hadrons and PF-bosons; 

• 1-loop electroweak radiative corrections and effective EW couplings by means of the DIZET 
v . 6 . 30 package; 

• higher-order leading-logarithm photonic corrections by means of the electron structure func- 
tions 123 H3 E3 EHj; 

• light pair corrections in the 0(a 2 L 2 ) leading-logarithm approximation including (optionally) 
the two-photon and singlet mechanisms. 

The code is applicable with the following restrictions: 

a) .©beam 3> me'- the energy has to be much larger than the electron mass; 

b) TOe/i/beam <C 9: extemely small angles are not described well, but the condition is fulfilled in 
practice for both small- and large-angle Bhabha measurements in the experiments at LEP, SLC 
and NLC; 

c) starting from the second order in a, real photon emission is integrated over, i.e. events with two 
photons separated from electrons are not generated. 

1 The codes are available upon request from the authors. 
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4.2 BHLUMI 



The Monte Carlo Program BHLUMI which has been used in the LEP analyses, is described in detail in 
ref. |22|. 

4.3 Comparison between BHLUMI and SAMBHA 

BHLUMI is compared with SAMBHA for integral and, for the first time, also differential distributions. 
The actual measurements are of calorimetric type. Therefore, event samples are generated with both 
programs, subjecting each event to a common set of calorimeter-like criteria (hereafter called CALO). 

In a first test the program codes were applied to the conditions established by the working group 
Event generators for Bhabha scattering |12| , with the result that all numbers were reproduced within 
the quoted accuracy. 

In a further test, about 10 8 Bhabha events were generated according to the calorimeter like con- 
ditions specified in sect. 5.1. This selection rejects a considerable part of events with real hard 
photon radiation. Therefore, the effect of mutual cancellation between virtual and real radiation is 
suppressed, which inevitably causes fairly large ^-dependent radiative corrections. The comparison is 
made quantitative in the form of the ratio 



p(t) 



^^sambha 

/At- der b hiumi/d£ 

d(7bhlumi/di 



and displayed in fig. [3] A linear logarithmic fit to the cross-section ratios and their statistical uncer- 
tainties gives 



p(t) = -(0.0039 ± 0.0002) - (0.0046 ± 0.0010) log- 



-t 



with (t) = -8.3GeV 2 . 



-t) 



log(-t) 

Figure 3: Cross section ratio p(t) as a function of log (—t), with t in units of GeV 2 . 

The two programs differ significantly, on average by 0.4%. At the present level of investigation, it 
cannot be excluded that there is a weak t dependence. 
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It is not so surprising to find a discrepancy for the differential quantity, while getting good agree- 
ment for the integral quantity. In fact, the key element is the far stronger restriction in the event 
selection for the two cases. In the integral case the events are accepted over the entire angular range 
of the luminometer, while for the differential analysis the same selection criteria are applied to a set of 
segments covering eventually the full range of the luminometer. This implies that an event accepted 
in the integral case is not necessarily accepted in the differential one owing to the more restrictive 
conditions, so that the sum of events accepted in the segmented luminometer is smaller than the 
number of events in the full luminometer. 



Table 2: Comparison between the codes NLLBHA (SAMBHA) and BHLUMI. Numbers are obtained by using 
conditions on table 19 in ref. ^21- The relative ratio 8r/r in per mille is defined by (YR-NOW)/YR. 
Last column gives the relative difference between BHLUMI (NOW) and SAMBHA(NOW). 



cut 


BHLUMI (YR) 


BHLUMI (NOW) 


5r/i 


NLLBHA (YR) 


NLLBHA (NOW) 


<5r/r 




0.1 


166.892 


166.879 


0.07 


166.948 


166.923 


0.14 


-0.26 


0.3 


165.374 


165.438 


-0.38 


165.448 


165.420 


0.16 


0.10 


0.5 


162.530 


162.616 


-0.52 


162.561 


162.25 


1.91 


2.25 


0.7 


155.668 


155.733 


-0.41 


155.607 


155.40 


1.33 


2.13 


0.9 


137.342 


137.425 


-0.60 


137.199 


137.32 


-0.88 


0.76 



For a quantitative understanding of this qualitative argument, a selection of events is presented 
as a function of the cut s ■ x c = s ■ x±X2, where xi is the fractional energy carried by the electron 
(or positron) (see table EJ. Obviously, a value for x c near to 1 selects configurations with small 
acollinearity, as opposed to cases with smaller x c , which favour larger acollinearity configurations. 
For a given opening angle, events with large acollinearity are hardly accepted; in other words the 
size of the cone opening angle defines the number of radiative events containing real emitted photons 
accepted or rejected. Consequently a larger or smaller final-state phase space is probed. Since virtual 
radiative contributions are unaffected by phase-space restrictions, the interplay between real and 
virtual radiative contributions strongly depends on the acceptance. The accuracy to which radiative 
corrections have to be treated becomes crucial. 

With the tight cuts required for the study of a differential quantity, as in the case investigated 
here, fine detailed aspects related to radiative contributions are necessarily probed. Therefore such 
studies open a new level of comparison between theory and experiment. 

5 Evaluation of the running in a simulated experiment 

Anticipating the application of the proposed method to measure the t dependence of a(t) on the 
data of a real experiment, a Monte Carlo simulation is carried out instead, in order to demonstrate 
the feasibility. An event sample is generated in the conditions of the DELPHI experiment using 
the existing program BHLUMI. In the next subsection the sample so obtained is confronted with the 
expectation of the new program SAMBHA. It should be noted that the t dependence of a(t), i.e. the 
quantity to be investigated, is stronger by about an order of magnitude than the possible differences 
in the intrinsic t dependences between BHLUMI and SAMBHA(see sect. 4. 3). 

5.1 Event generation 

The DELPHI detector and its performance are described in ref. - F° r the analysis, the relevant 
subdetector is the electromagnetic calorimeter STIC |30|. which covers the extreme forward and 
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backward directions. It has a ring structure with segmentation in both 9 and 4> covering \f—t ranges 
from 1.5 to 6 GeV for LEP1 energies and 3 to 12 GeV for LEP2 energies. 

Electrons, positrons and photons are observed as clusters. Their reconstruction is based on a 
cluster algorithm. The Bhabha events are characterized by two narrow high energy electromagnetic 
clusters opposite to each other and well inside the detector. The cluster algorithm is applied to the 
observed energy depositions in the cells of the electromagnetic calorimeter. Furthermore, the cluster 
with the highest energy satifies the more restrictive requirement to be at the radial position R between 
10 and 25 cm such as to cause no inefficiency for the opposite cluster. 

A Monte Carlo simulation has been performed using BHLUMI for three centre-of-mass energies of 
LEP: 91.2 (Z peak), 189 and 200 GeV. Assuming integrated luminosities / Cdt typical of the LEP 
experiments, the number of events passing the selection criteria is obtained and listed in table An 
event is attributed to ring z, if the highest energy cluster is reconstructed in this ring and the criteria 
listed below are satisfied. 

• Cluster reconstruction: 

The main criterion for merging adjacent cells is: 

V30mrady + V870 mrad ) < 
where the cluster centre is calculated as the energy-weighted cell centres. 

• Only the highest energy cluster in each hemisphere (referred to as F(jorward) and B (backward) ) 
is considered 

• Energy requirements: 
min(E F ,E B ) > 0.65 E hcmll 
max(E F ,E B ) > 0.94 E bcam 

This implies that the Bhabha events have not suffered from sizeable initial-state radiation effects. 

• Geometrical acceptance: 

The radial position R of the two opposite clusters must satisfy 

7cm < R f ,Rb < 28cm 

• Kinematics: 

The cluster centre and the nominal interaction point of the colliding e + e~ beams determine the 
dip angle 9. The quantity t is calculated from the dip angle 9 and the nominal centre-of-mass 
energy </s — 2E bcam according to 

t = ~s(l - cos6» max ) 

where # ma x is defined to be the dip angle of the cluster with the highest energy. 

The result of the Monte Carlo experiment is summarized in table [3] Ring 1 and ring 7 have been 
disregarded in order to exclude any inefficiency from border effects. 

5.2 Comparison and evaluation 

In this subsection the relevant observables and the parameters to be extracted are established and 
discussed. 

Each ring defines with its boundaries a bin (t m i n , t max ). The event numbers are to be equated 
to the corresponding theoretical prediction obtained from the formulae implemented in the program 
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Table 3: Numbers of events generated with BHLUMI 



\/i(GeV) 


91.2 


189 


200 


jCdt (pb- 1 ) 


75 


150 


200 


Ring 2 


1844850 


863571 


1028210 


Ring 3 


907754 


425586 


506131 


Ring 4 


513696 


240550 


286994 


Ring 5 


313218 


146731 


174740 


Ring 6 


201893 


94033 


112168 



SAMBHA. In order to extract the t dependence of a(t), eq. El is evaluated for each ring Ri defined by 
the geometry of the DELPHI luminometer. Equation then reads, for ring i : 



with the following definitions 



dt 

R > da 



°* = J "'dt 
a(U)\ 2 _ f R ' dt (a{t) xl 



a(Q) J J W - imin \ a (0) 



1 + An 



a(0)\ 2 a 



Table 01 contains the resulting theoretical values. 

Putting together the experimental and theoretical ingredients, i.e. the observed number of events 
Ni in each ring, together with the relevant luminosities J Cdt (from table^ and of , Ar^ (from table^J , 
we obtain the final formula: 

Ni 1 



a(0) J a? J Cdt 1 + An ' 
which can be exploited in a linear fit to access the parameters defining the t dependence of a : 

' ' ' ' ] ] - («o ± 5u ) + (til ± S Ul ) ■ log^r (9) 



M°)J (-*) 

The parameters of the fit are listed in table [SJ 



6 Discussion 

Table |S] demonstrates that for the case of the DELPHI setup (see sect. 5) and assuming typical 
integrated luminosities, the statistical accuracy is sufficient to verify the running of a for each of the 
three centre-of-mass energies. 

Equation [3] can be expanded in terms of Aa (see eq. |y. It is convenient to consider 

Ni 1 -U 

^TTA^ = no + nilog H> (10) 



n 



Table 4: Theoretical predictions for each ring of the three factors of eq. [7| For the conditions defined 
in sect. 5.1 the angular boundary of ring i is 6*i=arctan (7+3(i-l))/220). 



No. of ring 


1 


2 


3 


4 


5 


6 


7 










= 91.2 GeV 








2 


63.077 


C\ A TO O 

24.728 


12.170 


6.8694 


4.2517 


2.8120 


1.9552 


(a(ti)/a(oO 




1.0425 


1.0475 


1.0516 


1.0551 


1.0582 


1.0609 


1.0634 


1 + Ar z 




0.9426 


0.9440 


0.9412 


0.9395 


0.9240 


0.8915 


0.7982 










= 189 GeV 








2 


14.685 


5.7563 


2.8324 


1.5984 


0.9889 


0.6537 


0.4542 






1.0554 


1.0613 


1.0661 


1.0702 


1.0736 


1.0767 


1.0794 


l + An 




0.9377 


0.9390 


0.9360 


0.9329 


0.9165 


0.8858 


0.7898 










= 200 GeV 








2 


13.115 


5.1406 


2.5295 


1.4274 


0.8831 


0.5838 


0.4057 




1.0565 


1.0625 


1.0673 


1.0714 


1.0749 


1.0780 


1.0807 


l + An 




0.9376 


0.9387 


0.9352 


0.9330 


0.9158 


0.8847 


0.7896 










= 1000 GeV 








2 


0.5248 


0.2059 


0.1014 


0.0573 


0.0356 


0.0236 


0.0165 






1.0921 


1.0994 


1.1050 


1.1096 


1.1135 


1.1169 


1.1199 


i + An 




0.8622 


0.8620 


0.8590 


0.8545 


0.8398 


0.8084 


0.7205 










= 3000 GeV 








2 


0.0590 


0.0234 


0.0117 


0.0067 


0.0042 


0.0028 


0.0020 






1.1192 


1.1267 


1.1325 


1.1373 


1.1414 


1.1448 


1.1479 






0.8467 


0.8457 


0.8422 


0.8381 


0.8253 


0.7956 


0.6975 



rather than eq. [HI since in practice the integrated luminosity J Cdt is not known. The two coefficients 
no and n\ are obtained from a linear fit and contain the information on both the data and theory 
Their interpretation is : 

n = J Cdt ' fl + 2Aa((t» 

m = f Cdt - ( — -4 — r2Aa(t) 
J \dlog(-t) w 

The dependence on the integrated luminosity is given explicitly: obviously, one has n* = ttj ■ J Cdt by 
comparing eqs. |H1 191101 

In the ratio ni /no the dependence of the integrated luminosity drops out : 



d 



Aa = ^-(l + 2Aa((t)) 



dlog(-t) 2n V 

The slope dAa/dlog(— t), the quantity of interest, is then directly given by the ratio of the two 
experimentally measured quantities no and ni, namely ni/2no. The contribution of 2Aa((£)) is small 
with respect to 1 and can be neglected. The accuracy of the slope is determined by <5ni/2no, i.e. 
about 10% (see table EJ, which is far smaller than the absolute value of ni/2n . 
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Table 5: Table of fit results; the uncertainties Suq and 5u\ are uncorrelated. 





91.2 GeV 


189 GeV 


200 GeV 


U 


1.0573±0.0005 


1.0698±0.0008 


1.0703±0.0007 


Ul 


0.0242±0.0028 


0.0284±0.0041 


0.0318±0.0038 


("*) 


8.5 GeV 2 


36.6 GeV 2 


40.9 GeV 2 



On the other hand, no relates the integrated luminosity to Aa at the average value of t 

J l + 2Aa((t)) 

Making use of Aa((t)) as a priori knowledge the fitted no can be used to derive the integrated lumino- 
sity, which is the standard procedure. The statistical precision is given by Sno/no, which is of the 
order of 10 -3 . 

In addition, the hadronic contribution to Aa(t) (see fig-EJl may be deduced by subtracting the lep- 
tonic contribution, which is theoretically known precisely. The extraction of the hadronic contribution 
is only limited by the experimental precision. 

7 Conclusions 

A novel approach to access directly and to measure the running of a in the space-like region is proposed. 
It consists in analysing small-angle Bhabha scattering. Depending on the particular angular detector 
coverage and on the energy of the beams, it allows a sizeable range of the t variable to be covered. 

The feasibility of the method has been put in evidence by the use of a new tool, SAMBHA , to 
calculate the small-angle Bhabha differential cross section with a theoretical accuracy of better than 
0.1%. 

The information obtained in the t channel can be compared with the existing results of the s 
channel measurements. This represents a complementary approach, which is direct, transparent and 
based only on QED interactions and furthermore free of some of the drawbacks inherent in the s 
channel methods. 

The method outlined can be readily applied to the experiments at LEP and SLC. It can also be 
exploited by future e + e~ colliders as well as by existing lower energy machines. 

An extremely precise measurement of the QED running coupling Aa(t) for small values of t may 
be envisaged with a dedicated luminometer even at low machine energies. 
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